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ABSTRACT 

Dibaryons configurations are studied in the framework of the bound state 
soliton model. A generalized axially symmetric ansatz is used to determine 
the soliton background. We show that once the constraints imposed by the 
symmetries of the lowest energy torus configuration are satisfied all spurious 
states are removed from the dibaryon spectrum. In particular, we show that 
the lowest allowed state in the S = —2 channel carries the quantum numbers 
of the H particle. We find that, within our approximations, this particle 
is slightly bound in the model. We discuss, however, that vacuum effects 
neglected in the present calculation are very likely to unbind the H. 



1 Introduction 



Since it was first proposed by Jaffe [Q that some hexa-quark state could be stable against 
strong decays, great attention has been devoted, both theoretically and experimentally, to 
this issue. Jaffe's suggestion was based on a bag model calculation where he showed that 
color-magnetic interactions favor the existence of a stable flavour singlet S = —2 state 
(the so-called H-dibaryon). Later on, however, it was shown that the inclusion of effects 
neglected in Ref.|jl| (like i.e. symmetry breaking effects, center of mass corrections, pion 
cloud around the bag, etc.) tend to decrease the binding in a significant way 0, rendering 
it rather uncertain. Moreover, as recently discussed in Ref. 0, bag model predictions seem 
to be very sensitive to the bag constant which is not strongly constrained by empirical 
data. From the experimental point of view the situation is also unclear. Although some H 
weak decay events have been reported some time ago , other recent analysis provides no 
indication of a stable H-dibaryon New experiments are at the moment being carried 
out to investigate this issue further (see i.e. Ref. ). 

Theoretical calculations have also been performed using various other models, like i.e. 
lattice QCD, non-relativistic quark model and soliton models (for a rather extensive list 
of references, see Refs.[^ ^). Again, results have not been conclusive. In the case of the 
soliton models, most of the studies have been done using the collective coordinate SU{3) 
(see Refs.[^ ^ and references therein) extensions of the Skyrme model. In this paper we 
will use an alternative method based on the bound state approach^]. Since it has been 
recently shown |10|, |ll|] that within this scheme hyperon properties are remarkably well 
described, this will provide another interesting insight into the problem of the strange 
dibaryons stability. A previous attempt to investigate the structure of dibaryons within 
the bound state approach was done in Ref.lT^ where the H was found to be unbound. 
In that work, the simplified axially symmetric ansatz proposed in Ref. was used for 
the skyrmion field. Indeed, such a simplified ansatz predicts that the B = 2 soliton 
mass is twice larger than Msoi{B = 1) and is unstable. Here, we will use the improved 
axially symmetric ansatz proposed in Ref . |T^ . Although this ansatz does not correspond 
exactly to the lowest axially symmetric energy configuration, diskyrmion properties (like 
i.e. soliton mass, rotational energies, etc.) computed with it turn out to be very similar to 



those obtained with the lowest energy torus configuration numerically found in Ref. 1 15 



An intriguing point discussed in Ref.[0 is the existence in the bound state soliton 
model of states which are forbidden in the quark model. In fact, a similar situation 



was found in Refs.O, 14] for the case of non-strange dibaryons. In this paper we will 



show, however, that all these spurious states are removed from the spectrum once the 
constrains imposed by the symmetries of the problem are correctly taken into account. In 
particular, in our scheme the lowest S = —2 allowed dibaryon state is the flavour singlet, 
in agreement with the quark model prediction. 

This paper is organized as follows. In Sec. 2, we introduce the bound state model for 
arbitrary baryon number based on the improved axially symmetric ansatz. In Sec. 3, we 
show how to construct dibaryon wavefunctions which are consistent with all the symme- 
tries of the system. In Sec. 4, numerical results are presented and discussed. In Sec. 5, 
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conclusions are given. In Appendix A, we review the SU (2) sector of the model. Finally, 
in Appendix B the quantization rules and the dibaryon quantum numbers are discussed 
in detail. 



2 The Model 



We start with the effective action for the simple Skyrme model with an appropriate 
symmetry breaking term, expressed in terms of the SU (3)-valued chiral field U {x) as 



F2 

^Tr 
16 



d„Ud^'U^ 
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32e2 



Tr 



SB, 



where F^^ is the pion decay constant ( = 186 MeV empirically) and e is the so-called 
Skyrme parameter. In Eq.(|l|) F^^ is responsible for the explicit breaking of chiral sym- 
metry. We use the following form for F^^: 



SB 



48 

48 



Tr 



U + -2 



+ 



F^ml - F]^m\ 



24 



Tr 



(2) 



where A® is the eighth Gell-Mann matrix and and mx represent the pion and kaon 
masses respectively and Fk is the kaon decay constant ( = 1.22 F^). Eq.(^ accounts not 
only for the finite mass of the pseudoscalar mesons but also for the empirical difference 
between their decay constants. In previous calculations |^, 12 1 the kaon was found to 
be overbound to the soliton. It was recently shown |jTO| that this defect can be mostly 
eliminated if the difference in the decay constants is properly taken into account, as done 
in Eq.(l). Finally, Twz is the Wess-Zumino action. 



wz 



2407r2 



(3) 



which distinguishes between states with positive and negative strangeness. 

We proceed by introducing the Callan-Klebanov (CK) ansatz for the chiral field 



In this ansatz, Uk is the field that carries the strangeness. Its form is 



Uk = exp 

where K is the usual kaon isodoublet 

K 
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The other component, [/^, is the sohton background field. It is a direct extension to 
SU{3) of the SU{2) field u^, i.e., 



Ut, 

1 



(7) 



Replacing the ansatz Eq.(^ in the effective action Eq.(|I|) and expanding up to second 
order in the kaon fields we obtain the following lagrangian density for the kaon-soliton 
system 

C- = Csui2) + {D,K)\D^'K)-K^a^a^K 
2 



{k^K Tr {[a^.a'^f) - {D,Ky {D^K)TT{a,a'' 
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where Csu{2) is the effective pion lagrangian whose explicit expression is given in the 
Appendix A and 



is the baryon number current of the SU{2) configuration given by 

1 



B^' 



247r^ 



(9) 



(10) 



where 1^ = u^dyU-,,. 

In order to obtain the soliton background configuration we introduce the axially sym- 
metric ansatz 

= exp [ir ■ 7r„ F] (11) 

with 

Tin = sin B cos n(j) i + sin B surncj) j + cos B k (12) 

In Ref.[0 it was assumed that F = F{r) and Q = 6, where (r, 6, 0) are the usual spherical 
coordinates. As already mentioned, such an ansatz, which predicts 

R = Msoi {B = 2) / Msoi {B = 1) = 2.14 , (13) 

leads to an unstable soliton configuration. The lowest energy diskyrmion configuration was 
numerically found in Ref. ||15||. A very good variational approximation to such a solution 
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was proposed in Ref. |T^. In this case, F is a function of r only, but the variational 
function 

m 

Q = e + Y.9k sm{2ke) (14) 
fc=l 

is used for B. The coefficients Qk are determined by minimizing the soliton energy in 
the corresponding baryon sector. Using this ansatz one finds R = 1.94, which compares 
very well with the lowest energy solution value Rmin = 1-92. As already mentioned, other 
computed quantities like rotational energies, baryon radius, etc. are numerically also very 



close to those of Ref . . Since the use of the variational ansatz leads to a considerable 



simplification of the calculation, we will use it to describe the soliton background. The 



expressions corresponding to the SU{2) sector of the model have been obtained in Ref.|T4 . 
For completeness, they are summarized in Appendix A. 

Of course, for n = 1 the minimum soliton energy is obtained for = 0. Therefore, 
in this case the background soliton is still symmetric under combined spatial and isospin 
rotations A = / + / and the kaon field can be expanded in terms of spinor spherical 
harmonics, 

K{f,t) = k{r,t) yAiA,if) . (15) 
However, for n ^ 1 the background field is no longer invariant under A-rotations. In this 



case we use the consistent ansatz \12 



K{r,t) = k{r,t) f ■7r„(f) x , (16) 



where x is a two-component spinor. Note that, for n = 1, this ansatz reduces to Eq.(|T5D 
for the particular case A = 1/2, / = 1. These are precisely the quantum numbers of lowest 
kaon bound state when n = 1. 

Using the ansatze given above the explicit form of the kaon-soliton effective lagrangian 

is 

L = J drr^ {/(r) k^k - h{r) k^'k' + i\{r){k^k - k^k) - k'^k^rriK^ + K//)} , (17) 
where 



and 




Hr) = 1 + -^-^ ■ (19) 



The term linear in time derivatives, whose coefficient is 



A(r) = F' ^ , (20) 
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is due to the Wess-Zumino action, and 



eff 
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The diagonahzation of the hamiltonian obtained from the effective lagrangian Eq.(|T 
leads to the kaon eigenvalue equation 



--dr {r^hdr) +ml + Veff-fe^-2Xe 



k{r) = . 



(23) 



To obtain the hyperfine corrections to the dibaryons masses we proceed with the 
semiclassical collective coordinates quantization method, where the isospin and spatial 
rotations are treated as the zero modes [^. Then, we introduce the time-dependent 
spatial rotations R and the isospin rotations A such that 



K 



R Au^ A ^ 
RAK. 



(24) 
(25) 



The angular velocities in respect to the body fixed frame are given by 



(r~'r) 



ah 

A'^A 



^abc^c 
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(26) 
(27) 



Using ai and 02 as coefficients of the up and down spinor x in Eq.(|T6|) the substitution 



of Eqs.( 



in the full Lagrangian Eq.(||) yields 



-Msoi + Lk - T-Q + (Tiu;i + T2u;2)s2 + Tsc^s-si - {TiVti + T2VL2)t2 - T^Vt^^ti 



^X2(t^l + Lol) - Ii5n,l{VtiUJi + Vt2UJ2) + ]^I^{nVt^ - LO^f , (2^ 
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where Q 

si = 3 {2a4 di + [a^ — aia^)] ^2} , (29) 

S2 = 3|(l-a4)c?i + ^ [ai(l + "4) - as)] , (30) 

ti = n S2 ■, (31) 

t2 = 2(di + 4) 5n,i , (32) 



with di and (i2 given by 

di = 2£n I dr k*k 



-r^f cos^F/2 - _1^A(^2 ^'sinF) 
3 ' e^FK^dr ' 



2f r°° 4 
and the angular integrals 04 and 05 by 



(33) 

/ dr k*k- cos^ F/2 sin^ F , (34) 
JO 3 



1 /""^ 

q;4 = - / (i6' sin 9 sin^ 9 , (35) 
4 Jo 

as = r sin^ sin^ 6 6'^ . (36) 

JO 

In Eg. (|28|) is defined as T' = a*T-jaj and Xj are the moments of inertia of SU{2) sec- 
tor whose explicit expressions are given in Appendix A. The spin and isospin components 
Ji^ and /f'^ respectively, are calculated via 

J'f - — - ^ (37) 

Unlike the n = 1 case where the hedgehog symmetry of the skyrmion enforces the con- 
straint J + I = T, for n 7^ 1 only the constraint 

= -nil'/ + T3) (38) 

due to axial symmetry is left. 

The quantization of the rotational hamiltonian leads to 

+ l^yj + c,T,f + ^(4'T_ + fjT^) , (39) 

with the hyperfine constants Ci and C2 given by 

ci = 1 - (40) 

C2 = 1 - S2 . (41) 

Given this form of the rotational hamiltonian we should find the corresponding eigen- 
functions that satisfy the constraints imposed by the symmetries of the system. This is 
done in the following section. 

^In what follows we assume that Q = for n = 1. 
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3 Dibaryon wave functions 



In general the dibaryon wave functions are combinations of product states of the rota- 
tion matrices jb/(^) for isospin, for the angular momentum and the kaon 

eigenstates kT^{r-,t). Here Is and J3 are respectively the isospin and angular momentum 
projection on the lab. frame, I^^ and J^^ on the body fixed frame and T3 the projection 
of the kaon "spin" on the body fixed frame. Each of these product states have the form 



in) Dl^,,{u;) 



where 



-2(/; 



(42) 



(43) 



First we note that in these product states, for a given value of S, not all the values 



of isospin are allowed. This can be shown using the same method used in Ref . [p!6l for 
B = 1. Details are given in Appendix B. One obtains that for states in the minimal 
representations the allowed values of isospin are given by 



1 = ^ 



(44) 



where p < 6 + S should be odd if S is odd or even otherwise. This relation together with 
the axial symmetry constraint Eg. (^31) imply that J^^ is always even. 

To determine for each set of quantum numbers the correct linear combination of prod- 
uct states we should take into account all the symmetries of the problem. As discussed in 

2, the following 



Ref. |jT4| in addition to axial symmetry, the ansatz Eqs. 
reflection symmetries 



IJ) has, for n 



Ti2{-x,y,z) 



7r2[x, -y,z) 



TT2{x,y, -z) 




Ti2{x,y,z) 



(45) 
(46) 
(47) 



Denoting spatial rotations through an angle (j) about the z-axis by TZj (0) and isoro- 
tations through an angle 9 about the j-axis by T^jiO) the symmetry transformations are 
also generated by 



n({7r) 7^^(7^) 

V niin) 



(48) 
(49) 
(50) 
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The operator V denotes the parity operator which is defined as a space inversion and a 
sign change of all components of the pion field. 

The importance of the first two of these symmetry transformations on the construction 



of the dibaryon wavefunctions was first noticed in Ref. ||17||. In fact, as it discussed in detail 



in Sec.4-2c of Ref.|T^, for system with axial symmetry, the existence of an additional 
symmetry with respect to a vr-rotation about an axis perpendicular to the symmetry axis 
implies that this vr-rotation is part of the intrinsic degrees of freedom, and is therefore 
not to be included in the rotational degrees of freedom. We can express this constraint by 
requiring that the operator TZexti which performs a rotation TZ by acting on the collective 
orientation angles (external variables), is identical to the operator 7^j„i, which performs 
the same rotation by acting on the intrinsic variables. 



7^, 



ext 



int 



(51) 



Since in our case TZ can be either given by Eg. (^8|) or Eq.(^) we will have, in principle, 
two independent constraints. We will see, however, that when applied to the product 
states both symmetry operations produce the same result and therefore we are left with 
one single constraint. 

Then, in order to satisfy the constraints imposed by the symmetries of the ansatz we 
have to choose the following linear combination of product states: 



1 + (7^ex^)"'7^ 



int 



(52) 



where is a normalization factor. In writing Eg. (^2]) we have use the fact that acting 
on product states {TZextY = ("^mt)^- As we will see below this is satisfied by our product 
states. 

In order to have the explicit form of the dibaryon wave function we should apply the 
symmetry operators on the product wave function. For the collective operators we get 



7^f(7r)7^^(7r) 



-1 



-1 



(^) = (-) 



bf 



h~T. 



(53) 
(54) 



To determine the effect of the symmetry transformations on the intrinsic wave function 
we have to notice that apart from the contribution of the kaon field we have to include 
the contribution of the soliton itself. The latter contribution has been calculated in 
Refs. |1T7|, [T9[| . Using the fermionic nature of the B = 1 soliton configuration they foundQ 



7^^(7r)7^^(7^) ^^oi 



(55) 
(56) 



^Note that when one extends the hght flavour group from SU{2) to SU{3)^ as done in Ref.||T^ the 
minus phase is obtained from the Wess-Zumino term. 
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Finally, we have to calculate the effect of the symmetry transformations on the kaon field. 
Performing the symmetry operations on one kaon states we find 

n({TT)ni{n) kr-AS = ±1) = k^niS = ±1) , (57) 
7^^(7^)7^^(7^) kT,{S = ±1) = Ti k_T,{S = ±1) . (58) 



In deriving Eqs. (|57|j58|) , we have use the form Eq.(|T^) for 5 = +1 kaons and its charge 
conjugate for S* = — 1 kaons. Therefore, for a state with arbitrary value of strangeness S 
we obtain 

7^^(7^)7^^(7^) fcT3 = {-)-'^'k.T, , (59) 
7^f(vr)7^^(7r) A;t3 = {-)-'^'k.T, . (60) 

As we see, the only difference between the action of the symmetry transformation 

bf 



Eq.(48) and that of Eq . (|49|) is the presence of an extra phase {—)"^^' in Eq.(54). Since we 
have shown that in our case J^^ is always even, we see that, as stated above, both sym- 
metry operations produce that same effect when acting on our product states. Moreover, 



using the constraints Eqs. ( ^PD together with the effect of the symmetry operations on 



the intrinsic and collective wave functions Eqs.(p3-|60), it is easy to show that TZl^t = 



when applied to the product states. 

Therefore, the dibaryon wave functions that satisfy all the constraints imposed by the 
symmetries of the system have the structure 

\Ih,JJ,,S> = iV(Z}i,_2;,((^)Di^_^3(u;)fcT3(r,t)- 

(-)'^'-"/'^i3,2;.(^)^i-;.+T3(^)fc-T3(r,t)) , (61) 

where K = 1^ + T3. The normalization constant can be easily calculated. We obtain 



(2J + 1)(2/ + 1) 



N = . . . ^ ■ (62) 



2{1 + 5 jbf q5t3,o) 



57r 



To determine the parity of these wavefunctions we have to make use of the symmetry 



operation given in Eq.(^). In fact, as discussed in Sec.4-2f of Ref. ||18| the existence of 



such kind of symmetries implies that the parity operator V can be written as 

V = Sn-\ (63) 

where S (= VlZi^iji) in our case) acts on the intrinsic coordinates, while TZ^^ ( = [7^3(71)]^^ 
in our case ) acts on the collective variables. Using the fact that K is always integer and 
that 

S ipsoi = i'sol (64) 
S k^T, = i-f^' k^T, (65) 
T^ii^) Di,^2KiP)DiMK-T,){^) = (-)^(^-^^) Di,^2KmDl^iK-n)i^) (66) 
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we find that the parity of the wavefunction Eq. (|61|) is given by 



(67) 



From the explicit form of the wavefunction we note that for i^' = 0, T3 = only 
the "Fermi-allowed" combinations (satisfying the constraint (^—y+-^~^/'^ = —1 ) survive, 
whereas the "Fermi-forbidden" ones have zero norm. For the particular case of the two- 
nucleon system [S = 0), we have the constraint {—)^^'^ = —1 which is nothing but the 
generalized form of the Pauli principle. It is interesting to note that in our case this 
principle realizes in a more conventional way than in Ref.|T^ (see Eqs.(5.7a-b) of that 
reference). This is due to the more convenient choice of spatial coordinates we have used 
in the present work. In general, using the constraint Eq.(H) together with the explicit 
form of the dibaryon eigenstates Eq. (|6TD , one readily finds that all the spurious states 
obtained in Refs. |ll2|, 0, |l4l are removed from the spectrum. The quantum numbers of 
the allowed states with J < 2 are shown in Table 1. 

An important remark is that the wavefunction obtained above is not an eigenfunction 
of the rotational hamiltonian Hrot- Since the Coriolis term (last term in Eq. (|39|) ) does 
not commute with T3, the eigenf unctions of the Hrot will be combinations of those given 
in Eq.(^) with the same II3, JJ3, {J^^Y quantum numbers, but different values of T^. 



Using Eqs.(|| 
dibaryons we get 



3) one can finally obtain the dibaryon mass formula. For ground state 



(6^ 



where 5* is the strangeness of the state and Mrot{S) is the rotational contribution. Since 
for S" = and S = —1 only one value of is allowed (T| = for S" = and = 1/4 for 
S = —1), in these cases M^ot is simply given by the mean value of Hrot in the corresponding 
dibaryon state. For 5* = we get 



s=o 

rot 



1 

2T1 



j(j + i) - (j; 



+ 



1 

22o 



1(1 + 1) - 



+ 



(69) 



i.e., the same as 13 



For S = —1, the rotational contribution is given by 



M: 



s=-i 

rot 



+ 



1 

1 
2Jo 



Jij+i)-{j'/r 



/(/+l)-(40' + C2 




(40^ - V 



Cl 



+ 



i+j+ii 



(70) 



For S 



-2, two values of Tg^ ( = 0, 1) are in general allowed ^. Therefore, M^ot is 



given by the eigenvalues of 
1 

2Ji 



M. 



S=-2 
rot 



j(j+i) - if/r] + ^ 



/(/+i) + d 



+ 



21. 



Whether both values are aUowed for a given state depends on the values of the other quantum 
numbers. 
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I I 2X2 ' 2J3 X2V2 V ' ' ' '-'^ ' \ (^l) 



with M = ii' ± 1. In the construction of 2 x 2 matrix we have used the basis given in 
Eq.(|6TD ordered according to increasing T^. 

Note that although in some particular cases the expressions for the rotational correc- 



tions to S* = — 1 and S — 2 dibaryons agree with those given in Ref. |12|, in general they 



differ. This is due to the fact that in general the eigenfunctions used in that reference do 



not satisfy the constraints imposed by the symmetry transformations Eqs.( |i8| , ^9D . 



4 Results and discussion 

In our numerical calculations we will consider two sets of values for the parameters in the 
effective action. In one case (SET A), we consider the chiral limit in the SU{2) sector, 
m^r = and fit and e to reproduce the empirical and A masses. This corresponds 



to the result of Ref. 20 



F^ = 129 MeV, e = 5.45 . (72) 
The second set of parameters (SET B) is obtained for 7727^ = 138 MeV. It corresponds to 



the result of Ref. [21 



F^ = 108 MeV, e = AM . (73) 

In both cases we take the kaon mass and the ratio F^/F^^ to their empirical values 
rriK = 495 MeV, F^/F^^ = 1.22. When comparing our results with those of Ref. ||12||, it 



should be noticed that in that reference the ratio of decay constants was taken to be 1. 
As mentioned above this leads to an important overbinding in the B = 1 sector. The 
spectrum of non-strange and strange baryons obtained for our two parameters sets can 
be found for example in Ref.[]TU[ and it will be not repeated here. 



As shown in Ref. [|T^ when SET B is used and only gi is taken as variational parameter, 



the minimum in the B = 2 soliton mass is found for 

91 = -0.339 Qi^i = . (74) 

The inclusion of (72 as a second variational parameter does not lead to any significant 
improvement in the energy minimum. In addition we have also checked that for SET A 
the minimum appears almost at the same value of gi. For this reason we will use the set 
of gi given in (|7^ in all our calculations. 

In the Table 2, we present the calculated values of all the parameters appearing in the 
dibaryon mass formulae for both the massless and the massive pion cases. Non-strange 
sector parameters corresponding to SET B have already been given in Ref. . They are 
repeated here for completeness. The values of e indicate that kaons are less bound to the 
soliton than in the B = 1 case. Using the B = 1 kaon eigenenergies given in Ref. Jill, 
get 

As = e{B = 2) - e{B = 1) = 16 MeV (75) 
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for both, SET A and SET B. A similar result was found in Ref. ||12||, although in that case 
the value of Ae was somewhat smaller, Ae ~ 10 MeV. It should be noticed that when 
the empirical value of the meson decay constant ratio is used in order to eliminate the 
large kaon overbinding found in Ref.[jl2|, Ae becomes smaller. In fact, if we set all Qi = 
and use SET A as done in Ref.fT^ but keep Fk/F^^ = 1.22, we find Ae = 0. On the 
other hand, the use of the improved axially symmetric ansatz increases the value of Ae, 
the effect being larger for Fk/F^ = 1.22 (=16 MeV) than for Fk/F^ = 1 (= 8 MeV). 
Similar comments hold for SET B. 

The calculated rotational corrections M^ot to the dibaryon masses are shown in Table 
3. We list only the allowed states with M^ot < 250 MeV. For comparison, we also give the 
results of Ref. [0 and the sum of the rotational contribution to the lowest baryon-baryon 
state in each particular channel. We observe that our rotational energies are somewhat 
dependent on the choice of parameters, being smaller for SET B. For SET A they are in 
general very similar to those reported in Ref.|]l2|. An important exception is the lowest 
S = —2 state for which we predict roughly half of the value given there. The reason 
for this difference is mainly the smaller value of C2 obtained in our model. This is of 
some importance since in our calculation the rotational energy of this state lies below 
the corresponding threshold for both sets of parameters in contrast with the situation in 
Ref.|T^|. Another interesting S = —2 state to be discussed is the (1,1^) state. This is 



lowest state for which the off-diagonal terms in Eq.(|7lD do not vanish. Therefore, in order 
to calculate the corresponding rotational energy the 2x2 matrix has been diagonalized. 
It should be noticed than in doing so one goes beyond the 0{N~^) in the l/Nc expansion 
(where Nc is the number of colors)^. Strictly speaking this is not consistent with the fact 
that other contributions to the same order have been systematically neglected. However, 
if we just consider the corresponding diagonal term in Eq.(^) as the rotational correction 
to such an state one obtains M^oti^, 1"*") = 154 MeV (SET A). This rather large rotational 
correction is in disagreement with different calculations (see i.e. Ref. [§) that predict the 
existence of a low-lying state with these quantum numbers. 

Now we will focus our attention on the problem of the H stability. The different 
contributions to the corresponding binding energy are given in Table 4. As we see, for 
both sets of parameters, the H is bound in our model. Since the attraction found at 
the level of hyperfine (rotational) interactions is not enough to compensate the fact that 
kaons are less bound to the diskjnrmion, the main source of attraction turns out to be the 
rather low value of the B = 2 soliton mass obtained with the improved axially symmetric 
ansatz. A similar situation occurs in soliton models calculation based on the SU{3) 
collective coordinates approach. At this point it is important to recall that although the 
lowest energy state in the 5* = channel has the quantum numbers of the deuteron its 
identification with this particle it is not completely clear. In particular, a large binding 
( ~ ISOMeV ) is predicted for this state in the present model |T^. Since a considerable 
part of this overbinding comes from the 0{Nc) contribution to the dibaryon mass this 
might be an indication that the static axial symmetric torus does not provide the correct 

For a similar situation in a slightly different context see Ref.g 



12 



description of the B = 2 baryons. In using this configuration only the potential energy 
aspect of the dibaryon state is addressed, but all of the relative kinetic energy contribution 
has been neglected. In fact, it has been recently shown in Ref.[p^ where a bound two- 
skyrmion configuration was numerically studied on a discrete mesh, the torus state was 
formed only at the closest encounter of the two skyrmions, whereas most of the time the 
two-skyrmion system consisted out of two well separated baryons as expected for a "well- 
respecting" deuteron. Another effect that has been ignored in our calculation (and which 
is possible related with the lack of kinetic corrections mentioned above) is the Casimir 
effect due to zero point vibrations. A recent estimation of this effect using the SO (3) 
ansatz indicates that it tends to shift the mass of the B = 2 configuration relative to 
two B = 1 hedgehogs upwards. If we assume that this is the main source of repulsion that 
brings the deuteron mass to its physical value, then such a repulsion would be enough to 
unbind the H-dibaryon. 



5 Conclusions 

In this article we have studied the structure of strange dibaryons in the context of the 
bound state soliton model. In this approach, such dibaryons are assumed to be bound 
states of kaons and a B = 2 topological soliton. To describe the diskyrmion configuration 
we have used an axially symmetric ansatz where the dependence on the azimuthal angle 6 
is found through a variational method. Such an ansatz provides a very good approxima- 
tion to the numerically found lowest energy solution which also has axial symmetry. We 
have shown that once the constrains imposed by the symmetries of the torus background 
configuration are satisfied all spurious states are eliminated from the spectrum. In par- 
ticular, we obtained a generalized form of the Pauli principle for the lowest lying B = 2 
states with even values of strangeness. In the S = sector this implies that the lowest 
allowed state has the quantum numbers of the deuteron, while for the S = —2 sector the 
lowest state has the quantum numbers of the H particle. This is in agreement with the 
predictions of the quark based models. 

We have found that although kaons are less bound to the diskyrmion configuration 
than to a single soliton, the H-dibaryon is barely bound within our approximations. This 
binding is mainly due to the rather small diskyrmion mass with respect to two individual 
skyrmions. However, we know from the case of the deuteron (which is strongly overbound 
in the present model |]14[) that the static torus configuration tends to underestimate the 
B = 2 soliton mass. Moreover, numerical studies showed that the torus configuration 



is formed only at the closest encounter of two skyrmions. In this sense it would be 
interesting to see whether the dynamical departures from the lowest energy solution can 
be parametrized in terms of some collective coordinate. Zero-point fiuctuations of this 
coordinate could then provide a mechanism to increase the B = 2 soliton mass without 
affecting the B = 1 skyrmion mass. Another effect that is expected to decrease the 
binding energies (and which is probably very much related with the previous one) is the 
Casimir effect which leads to contributions of 0{N^) [^. Since these effects are expected 
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to give similar contributions to all the dibaryon states independently of their strangeness, 
it is clear that if they are strong enough to push the deuteron mass to its empirical value 
they are very likely to unbind the H. 
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Appendix A: The SU{2) Sector 

In this appendix, we give, for completeness the expressions for the SU{2) sector of our 



model. Most of them have been already published elsewhere ||14 
The lagrangian density Csu(2) in Eq.(H) is 

-hp^^m^^l-cosF) . (A.l) 

Replacing the modified ansatz given in Eq.([Tl|) the classical mass turns out to be a 
functional of F, 6 and $ which explicit expression is 



M[F, e, $] 



+ 



1 sin^F 

2g2 ^2 



f" + Q'' 



sin^ O sin^ F 



n/2 



sin'^ 



sin^e 2 ^,2 sin^F 
sm t7 



+ 



;i-cosF) 



with 



dF 
dr 



de 
'de 



and $' 



(A.2) 



(A.3) 



The minimization of Eq. (|A.2| ) leads to an Euler-Lagrange for each of the functions F, 
Q and $. The first, and more simple, is 



$" = . 



(A.4) 



Due to single- valueness of the chiral field we must have $ = mp, with n the baryon 
number. Using this result for the equation for B is 



n 



Ci sin e + C2— -r sin^ e 9" + Ci - C2n- 



'sin 9 



2\ 2 



sm 9 

— cos 9 = 0, 

smb* 



cos 9 9' 
(A.5) 



with 



C, = 2./ ir .n=F + and = 2. ^ ■ (A-^) 

As pointed out by Kurihara et al. [Q, if one uses the form 9 = ^ as done in Ref.|jl3| this 
equation implies 

Ci{l-n^)cose = . (A.7) 
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which can only be satisfied for n = ±1. Therefore, for n ^ ±1 there is a local unstability. 
Instead of solving Eq. ( [A.5|) numerically Kurihara et al. proposed to use the trial function 
defined in Eq.(p!4|). 

Finally, the equation for the chiral angle F is 



'^'^ 2 , -2 17// I -^1" 17/ 

r H — - sm i' i' + r 



—^di sin F cos F + ^F'^ sin F cos F 



- 2 



a;2 sin^ F cos F 



-r sin F = 



(A.8) 



with the boundary conditions F(0) = vr and F(oo) = 0. The explicit expressions of cti 
and a2 have been given in Eq . (|2^) . 

The use of the collective coordinate method for the quantization of the SU{2) La- 
grangian leads to the four last terms of Eq . (|28D . The explicit form of the moments of 
inertia Xj appearing in such an equation is 



dr sin^ F 



+ Xf" 



1 sin^ F 

C« + ~o lo Vi 



(A.9) 



where the coefficients Q and rji are given by 



,2sinf_e 
tan2 e 



= 7t£ de sine [Q" + n 
(2 = vr/o" rf^ sin^(l + cos2e) 
(3 = 2tt £ de sin e sin^ 6 
C4 = fvr , 



^1 
^2 

^4 



n^vr/o" sin^(l + cos2^)e'^ ''""'^ 
tt/o" c/^ sin^(e'^cos2e + 
27r/(f de sine B'^in^ 6 , 



2 sin^e 



'(A. 10) 



Appendix B: Quantization rules and the dibaryon isospin 

In case the background field is restricted to SU{2) one finds that in the bound state 
approach (as in the case of the normal SU{2) Skyrme model) there is a quantization 
ambiguity on whether a baryon number B configuration has to be quantized as a fermion 
or as a boson. Furthermore one has no information on which SU{3) multiplet a quantized 
bound state system with J, / and S should belong. In order to avoid these ambiguities 



in the bound state approach the authors of Ref. suggested to introduce a third light 
fiavour - degenerated with the u and d fiavours - called "funny strange" fiavour and to 
embed the bound state ansatz (4) into SU^A)^^'^ (where the label (/) distinguishes this 
group from a physical 5*^7 (4) fiavour group) as follows: 



Uif) = { yui'^ ^if) ^Jui^^ ) _ 



Now t/y) belongs to SU{3) 



if) 
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and U^J^ has in terms of K^-^^ the same form as Uk in terms of K, where 

Kif^ = ( = I ^0 I . (B.3) 




After collective quantization of the SU{3)^^^ flavour degrees of freedom (u, d and f) 
the Wess-Zumino term Q leads to a constraint on the "funny" right-hypercharge p!6| , 

yL" - ^ . (B.4) 

where A^^^ is the number of colours, B the baryon charge of the configuration and 5* the 
physical ( not "funny" ) strangeness. The derivation of this relation assumes only that the 
background field belongs to SU{2) and that K has the form (6). The result (|B.4|) is 
as well applicable for the S'?7(3)'''^^-extended axial rotor system with n > 2 where n = B. 
For B = 2 and iV^ = 3 we have therefore the constraint 



yi^) = 2 + I . (B.5) 

Of all possible 5[/(3)(^) multiplets with Y^^^ given by Eq.( [B.5| ) only the "minimal ones" 
(which are naturally the energetically lowest ones) will interest us. "Minimal" means 
that they are composed in a minimal way out of p SU{3)'^^^ triplet and q antitriplet 
representations. Therefore they have 

y^/) = f^ aud /j/' = |, (B.6) 

where /)j is the right isospin. 

The multiplets are then described by the SU{3Y-^^ rotation matrices 

^y?/)7(/)/(/)y(/)/(«/(« ^^l'^2,...,78) , (B.7) 

3 ^ R R R,z 

where 71,72,73 are of course the physical Euler angles for isospin rotations, Y'^^'^ is the 
funny hypercharge of a member of the multiplet, /^■^^ its isospin, Jg'^^ the projection of 
the isospin in the lab. frame and 1)^^ the projection of the right isospin in the body-fixed 
frame. 

Since the "funny strangeness" was introduced only to serve as an auxiliary quantity, the 
physical states should belong to those SU (2) submultiplets of a given 5f/(3)(-^) multiplet 
which do not have any "funny strange" component. In practice this means that we 
only work with the upper row of a given minimal S'f/(3)'''^^ multiplet. Thus we have 
Y = y'^-^'' = y^'^^ and / = Z^-^-* = /j/^ = p/2 where Y and / are the hypercharge and the 

^The derivation of this result is in complete analogy to the one for the quantization of an SU{3) 
skyrmion. In that case, however, the Wess-Zumino term leads to the constraint Y/j = A'c i?/3 0]. 
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isospin of the physical states. In this case, the 5*^7(3) '■■^^ rotation matrix ( |B.7| ) simphfies 

to 

which is equivalent, as the isospin content is concerned, to the usual SU{2) isospin rotation 
matrix 

'(71' 72, 73) (B.9) 

where is the isospin projection in the lab. frame and I^^ the one in the body-fixed 
frame. 

From Eqs.( |B.6] , p.5| ) we can now construct the allowed funny SU{?i)^^^ multiplets for 
a bound state configuration with B = 2 and given S. For S* = we find {p,q) = (0,3), 
(2,2), (4,1) and (6,0) which correspond to the usual 6 quark flavour multiplets: 10, 27, 
35 and 28. For S = —1, we find {p,q) = (1,2), (3, 1) and (5,0) which correspond to the 
SU{3)^^'^ multiplets: 15, 24 and 21 respectively. Note that these numbers are the same 
as for the 5 quark SU{3) multiplets which follow after removal of one quark from the 6 
quark system. 

In Table ^ all minimal "funny" multiplets {p, q) for the B = 2 bound state configura- 
tions with strangeness S between and —6 are listed. From that table one can deduce 
that the states characterized by the rotation matrix of the background field (|B.8| ) corre- 
spond in the quark model to those minimal quark flavour configurations which involve 
only the u and d quarks of a given B = 2 state. E.g. for S = —6 the funny multiplet 
is a singlet signalling that there is no u, d content. The role of the strange quarks (6 in 
the S = —6 system) is taken over by the kaons coupled to the n = 2 soliton background. 
Note that any B = 2 state with S < —7 or S* > 1 corresponds to a non- minimal SU{3)^^^ 
multiplet for the rotor. 
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Table captions 



Table 1: The allowed quantum numbers in the B = 2 bound statesystem. In the rows 
the following quantities are listed: (1) the strangeness S, (2) the allowed values of isospin 
/, (3) the magnitude of the isospin projection in the body fixed frame, < /, (4) the 
magnitude of the projection of the kaon grand-spin on the body fixed frame, IT3I < |>S'/2|, 
(5) the magnitude of the projection of the total angular momentum on the body fixed 
frame, | Jg^l = 2I/3 + T3I (see Eq.(|H|)) with \J^^\ < 2, the parity P as given by Eq.(§7p, 
and the number of states with the same J^^ quantum number, (6) the total angular 
momentum J < 2 and the corresponding SU{3) representations. 



Table 2: Numerical values of the parameters appearing in the dibaryon mass formulae 
Eqs.( pE| - |7TD . SET A corresponds to the massless pion case while SET B corresponds to 



massive pions. 



Table 3: Rotational contributions to the dibaryon masses. Listed are those corresponding 
to the allowed states with Mrot < 250 MeV. SET A and B are as in Table 2. A^A^, A^A 
and AA stand for the sum of the rotational contributions to the corresponding particles 
and serve as rotational threshold in each particular channel. 



Table 4: Contributions (in MeV) to the mass of the H-particle given relative to twice 
the corresponding contribution to the A mass. 



Table 5: Allowed "funny" multiplets {p,q) for the B = 2 bound state configurations. 
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